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Abstract 

We prove six theorems concerning exponentially accurate semiclassical quantum 
mechanics. Two of these theorems are known results, but have new proofs. Under 
appropriate hypotheses, they conclude that the exact and approximate dynamics of 
an initially localized wave packet agree up to exponentially small errors in h for finite 
times and for Ehrenfest times. Two other theorems state that for such times the wave 
packets are localized near a classical orbit up to exponentially small errors. The fifth 
theorem deals with infinite times and states an exponentially accurate scattering result. 
The sixth theorem provides extensions of the other five by allowing more general initial 
conditions. 
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1 Introduction 



This paper is devoted to proving several theorems concerning exponentially accurate approx- 
imations to solutions of the time-dependent Schrodinger equation 

d 

ih — 'i!{x,t,h) = - — A<iJ{x,t,h) + V{x)'^{x,t,h) (1.1) 

in the semiclassical limit — > 0. 

The semiclassical approximation of quantum dynamics has been the object of several 
recent investigations from different points of view. One approach uses coherent state initial 
conditions and approximates the evolved wave packet by suitable linear combinations of 
coherent states. Another approach considers the Heisenberg evolution of suitable bounded 
observables and approximates the corresponding operators by means of Egoroff 's theorem. 
The goal of both approaches is to produce accurate, computable approximations as h goes 
to zero, for as long a time interval as possible. In scattering situations the time interval is 
the whole real line. 

There are several results concerning the propagation of certain coherent states for finite 
time intervals. Early results fl^, |^ constructed approximate solutions that were accurate 
up to 0{h^/^) errors. Later approximations were constructed with 0(^''^^) errors for any / 
P, H, [1^ . Very recently, approximations were constructed with errors of exponential order 
0{e-^/^) with r > in Inl (see also 



The validity of the corresponding approximations for time intervals of length 0(ln(l/^)), 
the so-called Ehrenfest time-scale, has been established up to 0{h}^'^) errors in and up 
to 0{e^/^°') errors with < a < 1 in |]11| . There is physical intuition and evidence that the 



Ehrenfest time scale is the natural limit for the validity of coherent state type approximations. 
This issue is studied in detail for the quantized Baker and Cat maps in 0. 

Approximations have been constructed for infinite times in the context of scattering 
theory for coherent states. Approximate solutions with errors of order 0{h^^'^), uniformly in 
time, are produced in [0. This yields approximations for the scattering matrix with errors 
that are also 0{h^^'^). Related results for another class of states can be found in ^ . 



Corresponding results for the approximation of observables in the Heisenberg picture can 
be found in ||T^ for approximations with Oifi^^"^) errors for any I for finite times. Approx- 
imations with exponentially small errors both for finite times and for Ehrenfest times are 
constructed in fl]] and 0. 

The exponentially accurate results mentioned above, and those we present below, are 
obtained for Hamiltonians that satisfy certain analyticity conditions. The approximations 
are generated by optimal truncation of asymptotic series. 



For further information, we refer the reader to the review articles IS, 15, 12 . 



The present paper is concerned with the propagation of coherent states in the spirit of 
the first approach described above. We present a new construction of approximate solutions 
to the time dependent Schrodinger equation that is an alternative to the one presented in 
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The new expansion has several advantages. In addition to being exponentially accurate 



up to the Ehrenfest time scale (Theorems jST^ and p.3| ), it allows us to extend our previous 
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results in four separate directions: 

1. We get exponentially precise localization properties for both the approximation and the 



exact solution for both finite times and Ehrenfest times (Theorems 3^ and 

2. We get exponentially accurate information on the semiclassical limit of the scattering 
matrix for suitable short range potentials (Theorem |3.5| ). 

3. The new algorithm is superior for numerical computation. The work done to construct 
the approximate wave function for one value of h is used for the construction for all smaller 



values of h. This should be contrasted with the construction of |11] where every calculation 
must be redone for each value of h. 

4. The results in |]ll| concern the propagation of initial coherent states given by a linear 
combination of a finite number of elementary coherent states. We can control the new 
approximation as a function of A^, which also allows us to extend the validity of all previous 

In this case however, 



results to a more general set of initial states (Theorem 3^ 



the 



algorithm requires the computation of different quantities as h varies. 



The technical difference between the present construction and the one in [TJ] is the 
following: In both papers, we use a suitable time dependent basis to convert the PDE 
( p..l| ) into an infinite system of ODE's for the expansion coefficients in that basis of the 



solution to the Schrodinger equation. In |Tl|] we then construct the approximate solution by 
approximating this infinite system by a finite system, which we solve exactly. In the new 
approach, we substitute an a priori expansion in powers of ^^/^ into the original infinite 
system of ODE's. We construct our approximate solution by keeping a finite number of 
terms. This turns out to be quite efficient. 

The new approximation also plays a vital role in the construction of an exponentially 



accurate time-dependent Born-Oppenheimer approximation 13 . 



2 Coherent States and Classical Dynamics 

We begin this section by recalling the definition of the coherent states (pjiA, B, h, a, t], x) 

described in detail in ||10|. A more explicit, but more complicated definition is given in 0. 

We adopt the standard multi-index notation. A multi-index j = (ji, j2, . . . , jd) is a 

d-tuple of non- negative integers. We define \j\ = Y.t=i jk, = ■ ■ - x-^f, 

alii 



j! = (Ji!)(j2!) ■ ■ ■ (jd!), and = ^^^^^^^ . 

Throughout the paper we assume a e IR'^, rj e IR'^ and h > 0. We also assume that A 
and B are d x d complex invertible matrices that satisfy 

A^B - B^A = 0, 

A* B + B* A = 21. (2.1) 

These conditions guarantee that both the real and imaginary parts of BA^^ are symmet- 
ric. Furthermore, IieBA~^ is strictly positive definite, and (ReBA^^) ^ = A A*. 

Our definition of ^Pj{A, B, h, a, rj, x) is based on the following raising operators that are 
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defined for m = 1, 2, . . . , 

Am{A,B,h,a,riy = 



2h 



d d Q 

n=l n=l ^^-n 



Tlie corresponding lowering operators Am{A, B, h, a, rj) are tlieir formal adjoints. 

These operators satisfy commutation relations that lead to the properties of the 
(pj{A, B, h, a, T], x) that we list below. The raising operators Am{A, B,h,a,r])* for 
m = 1, 2, . . . , d commute with one another, and the lowering operators Am{A, B,h,a,ri) 
commute with one another. However, 

Am{A,B,h,a,r])* An{A,B,h,a,r]) - An{A, B,h,a,r]) Am{A, B,h,a,T]y = 6m,n- 

Definition For the multi-index j = 0, we define the normalized complex Gaussian wave 
packet (modulo the sign of a square root) by 

MA B, h, a, r/, x) = 7r-'^/^r^/^(det(v4))-i/2 

X exp { - ( (a; - a), B A'^ {x - a) ) / {2h) + i [r], {x - a)) /h} . 
Then, for any non-zero multi- index j, we define 

MA,B,h,a,rj,-) = -^ iAM,B,h,a,v)T (A(A5,^,a,r/)*r^ ••• 

X {AdiA,B,h,a,r]yy'' , B, h, a, t], ■). 

Properties 1. For A = B = I,h = l, and a = 7] = 0, the (pjiA, B, h, a, 1], ■) are just 
the standard Harmonic oscillator eigenstates with energies \j \ + d/2. 

2. For each admissible A, B, h, a, and t], the set { 4>j{A, B, h, a, t], ■) } is an orthonormal 
basis for L'^{IR'^). 



3. The raising operators can also be given by another formula that was omitted from ||10 
in the multi-dimensional case. If we set 

g{A, B, h, a, x) = exp | — ( (x — a), (^BA~^^ {x — a))/ {2%) — i{rj, (x — a) )/^ | , 

then we have 

(^„(A,5,^,a,r/)>) (x) 

1 d _ g 

^nm ^ {g{A, B, h, a, x) tpix) ) . 



2 g{A, B, h, a, x) dxn 
4. In [0, the state (f)j{A, B, h, a, rj, x) is defined as a normalization factor times 

HjiA; h'^/^ \A\^^ {x - a)) ^A, B, h, a, r/, x). 

Here 7ij(v4; y) is a recursively defined IjI**^ order polynomial in y that depends on A only 
through Ua, where A = |y4| f/4 is the polar decomposition of A. 
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5. By scaling out the |y4| and h dependence and using Remark 3 above, one can show that 
'Hj{A] y) e~y is an (unnormahzed) eigenstate of the usual Harmonic oscillator with energy 
\j\ + d/2. 

6. When the dimension d is 1, the position and momentum uncertainties of the 

(f)j{A, B, h, a, T], ■) are J {j + 1/2)% \ A\ and J{i + 1/2)% \B\, respectively. In higher dimen- 



sions, they are bounded by y (|j| + d/2)% and ydil + d/2)% \\B\\, respectively. 

7. When we approximately solve the Schrodinger equation, the choice of the sign of the 
square root in the definition of 0o(^) -B, %, a, i], ■) is determined by continuity in t after an 
arbitrary initial choice. 

8. We prove below that the matrix elements of [x — a)™ satisfy 

I ( 4>j{A, B, %, a, T], x), (x — a)™ B, %, a, r], x)) \ 



and 



< ;iH/2 (72d)l™l p||H ^{\k\ + l){\k\ + 2) • • ■ (|A;| + |m|). 



{4>j{A, B, %, a, r], x), {x - a)"^ (j)k{A, B, %, a, r], x)) = 0, if 



> \m\. 



(2.2) 



We now assume that the potential V : IR'^ — >■ iR is smooth and bounded below. Our 
semiclassical approximations depend on solutions to the following classical equations of mo- 
tion 



(2.3) 



d(t) 


= vit), 




= -VViait)), 


Ait) 


= I Bit), 


Bit) 


= tV^^\ait))Ait), 


Sit) 


= "f - V(a(t)) 



where V^"^^ denotes the Hessian matrix for V, and the initial conditions A(O), 5(0), a(0), 
r/(0), and 5(0) = satisfy i^. 

The matrices Ait) and Bit) are related to the linearization of the classical flow through 
the following identities: 



Ait) 
Bit) 



^'^'^ AiO)+^^BiO), 



daiO) 
dr]it) 



BiO) 



driiO) 
. dr]it) 



AiO). 



d7]i0) ' ' 9a(0) 

Because V is smooth and bounded below, there exist global solutions to the first two 
equations of the system ( |2.3| ) for any initial condition. From this, it follows immediately 
that the remaining three equations of the system (p.3|) have global solutions. Furthermore, 
it is not difficult ||, 0] to prove that conditions (|2.1|) are preserved by the flow. 
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The usefulness of our wave packets stems from the following important property [0. If 
we decompose the potential as 

V{X) = Wa{x) + Z,{X) = Wa{x) + {V {x) - Wa{x)), (2.4) 

where Wa{x) denotes the second order Taylor approximation (with the obvious abuse of 
notation) 

Wa{x) = V{a) + 1^(1) (a) (x — a) + V^'^\a) {x — a)^/2. 
then for all multi- indices j, 
d 



e^"^'^^ <P,{A{t),B{t),h,a{t),v{t),x) 

A + Wa^t^ix)] [e^^W/'^</.,(A(t),5(t),^,a(t),77(t),a;)l , (2.5) 



if A{t), B(t), ait), rjit), and S(t) satisfy (|2.3|). In other words, our semiclassical wave packets 
ipj exactly take into account the kinetic energy and quadratic part Wa(t) (x) of the potential 
when propagated by means of the classical flow and its linearization around the classical 
trajectory selected by the initial conditions. 

In the rest of the paper, whenever we write (j)j{A{t), B(t),h,a(t),ri(t),x), we tacitly as- 
sume that A(t), B{t),a(t),ri{t), and S(t) are solutions to (|2.3|) with initial conditions satis- 
fying (P). 



3 The Main Results 

In this section, we list our results concerning the propagation of semiclassical wave packets. 
The first is the construction of an approximate wave function that agrees with the exact wave 
function up to an exponentially small error. The construction is quite explicit. It depends 
on the somewhat arbitrary choice of a parameter g > 0. 

The precise result is summarized in the following theorem: 



Theorem 3.1 Suppose V{x) is real and bounded below for x E IR^ . Assume V extends to 
an analytic function in a neighborhood of the region = {z : \ Im Zj\ < 6} and satisfies 
\V{z)\ < M exp (t|2;P) for z E Ss and some positive constants M and r. 

Fix T, choose a classical orbit a{t) for < t < T , and consider an arbitrary normalized 
coherent state of the form 

ij{x,0,n)= J2 c,(O)0,(A(O),5(O),/i,a(O),r/(O),x). 

\3\<J 

There exists a number G > 0, such that for each choice of the parameter g G (0, G), there 
exists an exact solution to the Schrddinger equation, 

ih— = - — A^f + V"^, 
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with \l/(x,0,^) = tlj{x,0,h), that agrees with the approximate solution 

ij{x,t,h) = e^^W/^ c,{t,h)(^,{A{t),B{t),h,a{t),r]{t),x), 

\j\<J+3g/h-3 

up to an error whose L^{IR'^) norm is bounded by C exp { —'jg/h }, with 7^ > 0. Furthermore, 
the complex coefficients Cj{t, Ti) are determined by an explicit procedure. 

The second result shows that the approximate wave function of Theorem is concen- 
trated within an arbitrarily small distance of the classical path up to an exponentially small 
error if g is chosen sufficently small. 



Theorem 3.2 Suppose that the hypotheses of Theorem are satisfied and that b > is 
given. For sufficiently small values of the parameter g > 0, the wave packet ip{x,t,h) is 
localized within a distance b of a{t), up to an error exp { —Tg/h}, with Tg > 0, in the sense 
that 



'\x-a{t)\>b 



{x,t,h)\ dx 



< exp {-Vg/h} . 



Next, we turn to the validity of the approximation and its localization properties on the 
Ehrenfest time scale, i.e. when T is allowed to increase with h as \\i{l/h). 



Theorem 3.3 Suppose the assumptions of Theorem \3. 1\ are satisfied except that the upper 
bound on V is replaced by \V{z)\ < M exp {t\z\) for z & Ss and some positive constants M 
and T. Further, assume the existence of a constant N > and a positive Lyapunov exponent 
A so that \\A{t)\\ < N exp{Xt), for all t > 0. Then, for sufficiently small T' > 0, there exist 
constants C > 0,7' > 0, a E (0, 1), and a' G (0, 1), and an exact solution to the Schrddinger 
equation that agrees with the approximation 

ij{x,t,h) = e^^W/" Yl c,{t,h)(P,{A{t),B{t),h,a{t),r]{t),x), 
\j\<J+3/n''' -3 

up to an error whose norm is bounded by C exp { /h'^ }, whenever <t < T'ln{l/h). 
Moreover, ifr can be taken arbitrarily small, we can chose T' = -^(l—e) where e is arbitrarily 
small. 



Remark: The semiclassical approximation of observables in the Heisenberg picture holds 
for any T' < 2/(3A), when r << 1, as shown recently in That time interval is longer 
than those for which a localized coherent state can approximate the evolution of an initial 
coherent state, which is characterized by T' < 1/(2A). See for a study of related issues 
on quantized hyperbolic maps on the torus. 
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Theorem 3.4 Suppose the hypotheses of Theorem are satisfied and that b > is given. 
Then, for sufficiently small T' > 0, there exist T' > 0, a E (0, 1), and a' G (0, 1), such that 
the approximation of Theorem pTl satisfies 



I \il){x,t,h)\^dx\ <exp{-r7^''} 

J\x~a{t)\>b I 



whenever < t < T' 

Moreover, ifr can be taken arbitrarily small, we can chose T' = ^(1— e) where e is arbitrarily 
small. 

We also explore the validity of the approximation in a scattering framework and its con- 
sequences on the corresponding semiclassical approximation of the scattering matrix S{h). 
This requires assumptions on the decay of the potential and its derivatives at infinity. 

For scattering theory, we assume V satisfies the following decay hypothesis. 

D: There exist (3 > 1, Vq > 0, and Vi > 0, such that for all x G IR^ and all multi-indices 
m G 

D-Vi.) < (3.1) 



where {x) = ^/l + 



Theorem 1.2 of shows that under the hypothesis D, the solution of the classical equa- 
tions ( |2.3| ) satisfies the following asymptotic estimates: 

For any a_ G iR^ ^ G IR'^ such that (a_, r/_) G IR^'^\S, where S C {(a„, r/„) G 
IR^'^ : 7^ 0} is closed and of Lebesgue measure zero in IR^, there exists (a+, rj^) G IR^'^, 

7^ 0, and 5*+ G iR such that 

lim \a{t) — a± — ri±t\ = 0, 

t — ^ioo 

lim \T]{t) - ?7±| = 0, 

I— >±oo 

lim \S{t)~tr]^_/2\ = 0, 

t—^ — 00 

lim \S(t)- S+-tr]l/2\ = 0. (3.2) 



Moreover, for any d x d matrices {A_, B^) satisfying condition ( p. ID , there exist matrices 
(A+, B+) G Md{(rf satisfying (|A|), such that 

lim \\A{t) - A± -iB±t\\ = 0, 

t—*±oo 

lim \\B{t)-B±\\ = 0. (3.3) 
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Our assumption D implies that V is short range. It follows that if Ho{h) = — 2~^' 
then the wave operators defined by 

n'^{h) = s- lim emn)s/n ^-rHo{h)s/h (^34^, 

exist and have identical ranges equal to the absolutely continuous subspace of Hifi). As a 
result, the scattering matrix 

S{h) = Q~{hyQ+{h) (3.5) 

is unitary. 



Theorem 3.5 Suppose d > 3 and assume hypothesis D. Let (a_, ?]_) G IFi?'^\£ and 
(A_, B_) G Md{(I!f satisfy condition (^). Let Cj(-oo) G (U, for j G N'^, with \j\ < J, 
such that J2\j\<j |cj(— cxd)P = 1. Then, there exist (a+, rj^) G IR^'^, 5+) G Md{(DY sat- 
isfying (\2. 1\ ), 5*+ G iR and explicit coefficients Cj{+oo, h) G (C , for all j G N'^, with \j\ < Jn 
with Jn = J + 3g/h — 3 such that for some 'y>0, C>0, g>0 (depending on the the 
classical data), the states defined by 



^^{A^, B^,h,a_,ri_,x) = Cj{—oo)(t>j{A_,B^,h,a_,ri_,x) 

\j\<J 

<^+iA+,B+,h,a+,r]+,x) = e*^+/^ ^ Cj{+oo,h)(Pj{A+, B+,h, a+,r]+, x) (3.6) 

\j\<Jn 

satisfy 

||5(;i)<l>_(A_,5_,^,a_,r/_,-) - $+(A+,5+,;i,a+,r/+,-)|L2(K^) < C e-^/^ 
if h is small enough. 

Finally, we address the question of the generalization of the initial coherent state, whose 
evolution can be controlled up to exponential accuracy in the different settings considered 
above. 

For {a,ri) G IR^'^, we define Kh{a,r]) to be the operator 

(A,(a,r7)/)(x) = h-'l^e^^^^^^-^^^/'^ f{{x-a)/Vn). 
We define a dense set C in L'^{IR'^), that is contained in the set S of Schwartz functions, by 

C = |/(x)=^c,0,(Z^i^ 1, 0, 0, x) G 5, 

such that > with ^ |cjf < e"^*', for large J \ . (3.7) 

lil>^ J 
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Remark It is easy to check that the inequahty in ( |3.7| ) is equivalent to the requirement 
that the coefficients of / satisfy 

for large Another equivalent definition of C is 

C = U>oe-*^''°5, 

where = — A/2 + is the harmonic oscillator Hamiltonian. The set C is sometimes 
called the set of analytic vectors for the harmonic oscillator Hamiltonian. 



Theorem 3.6 All theorems above remain true if the initial condition has the form 

ilj{x,0,h) = {Ania,ri)(p){x), 

where (p E C. 



Theorem |3.1| is proved in Section ^. Theorem p.2| is proved in Section |^. Theorems p.3| 
and ^]4| are proved in Section ^. Theorem |3.5| is proved in Section ^ Theorem |3.6| is proved 



in Section [I^. 

4 An Alternative Semiclassical Expansion 

1 /2 

In this section we derive an expansion in powers oih ' .In later sections we perform optimal 
truncation of this expansion to obtain exponentially accurate approximations. 

We wish approximately to solve the equation 

dib 

th-£ = -yA^ + Vix)^lj, (4.1) 

with initial conditions of the form 

tlj{x,0,n) = ^ Co,,(O)0,(A(O),5(O),/i,a(O),r/(O),x), (4.2) 
\j\<J 

where ^ |coj(0)|^ = 1. 
\j\<J 

We can write the exact solution to this equation in the basis of semiclassical wave packets, 

ij{x,t,n) = e'^^'^/^ (^At,h)(f)j{A{t),B{t),h,a{t),7]{t),x). (4.3) 

j 

Note that the sum is over multi- indices j. The infinite vector c whose entries are the coeffi- 
cients Cj satisfies 

ihc = K{t,n)c, (4.4) 
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where K{t,h) is an infinite self-adjoint matrix. 

The matrix K{t,h) has an asymptotic expansion in powers of h^^'^. The cubic term in 
the expansion of V{x) around x = a{t) gives the leading non-zero term of order h'^^'^. The 
quartic term in the expansion of V{x) gives the term of order h'^^'^, etc. Thus, we can write 

oo 

K{t,h) ~ ^'^'K,{t), (4.5) 

fc=3 

with 

K,it) . E (iql^ XW". ,4.6) 

|m|=fc 



where X[t)"^ is the infinite matrix that represents (x — a)"^. Explicit formulas \TU\ 

show that entries of X{t)"^ and Kk{t) do not depend on h. 
We formally expand the vector c as 

cit,h) = Co{t) + h^'^ci{t) + ^^c^it) + ... 

= (4.7) 

k 

We denote the j**^ entry of Cfc(t) by Ckj{t). Note that k is a non-negative integer, and j is 
a multi-index. We substitute the two expansions ( [4.5| ) and ( [4.7| ) into ( [4.4| ) and divide by h. 
We then equate terms of the same orders on the two sides of the resulting equation. 

Order 0. The zeroth order terms simply require 

ico = 0. (4.8) 
From (|i^), the solution is obviously Coj(t) = Coj(O). We note that Cojit) = if |j| > J. 

Order 1. The first order terms require 

zci = Ks{t)co{t). (4.9) 

We solve this by integrating. Because of ( |4.2|) , Ci(0) = 0. From the form of Co(t), only 
finitely many of the entries of c-i(t) are non-zero, and Ci o = whenever |?'| > J + 3. In c? 

.pace.— ..WKasa..o.(^--^)_e..es. 

Order 2. The second order terms require 

= K^it)coit) + K,{t)ciit). (4.10) 
Again, we can solve this by integrating with C2(0) = 0. The only entries of C2{t) that can 
be non-zero are C2j{t) with \j\ < J + Q- In d dimensions, there are at most ^ ~^ ^ ^ 
non-zero entries. 



11 



Order n. In general, the rfi^ order terms require 

n-l 

iCn = J2 ^n+2-k{t)Ck{t). (4.11) 
fc=0 



To solve this, we simply integrate. We observe that Cnj{t) can be non-zero only if |j| < J+3n. 

3n 
d 



In d dimensions, there are at most ( "I" "I" i non-zero entries. 



Our expansion is different from the one constructed in 11 , and it is different from the 



Dyson expansion used in |[TT|]. All three of these expansions are asymptotic to the exact 
solution of the Schrodinger equation. We note that the main construction in yields 
a normalized wave function. The expansion derived above does not generate normalized 
wavepackets. 



To prove that this expansion is asymptotic, we apply Lemma 2.8 of [|T0|. To check the 
hypotheses of that lemma, we do the expansion above through order (Z — 1) to obtain Co(t), 
Ci{t), . . . , We substitute these into (|4.7|) with the sum cut off after k = I — 1. We 

then use the result in ( ^73|) and compute 

^i{x,t,h) = th-^{x,t,h) + —A^{x,t,h) - V{x)^{x,t,h) (4.12) 



Because the Cfc(t) solve ([4. 81), ( [4.9| ), ( [4.10| ), etc., there are many cancellations. We obtain 

^==0 |j|<j(o 



(4.13) 

Here, J{1) = J + 3/ — 3, and for each g, wj^'^j-^{x) denotes the Taylor series error 
W^iS(^) = V{x)- E ^^^^f^^ (X - aityr 

|m,|<g 

\m\=q+l 

for some Cmix,a{t)) = a{t) + 6'^,^.,a(j) (x - a{t)), with ft^^^^^w e (0, 1). 

If V is C'-^'^ on some neighborhood of {a{t) : t e [0, T] }, then each w'f(l^{x) that 
occurs in (|4. 13| ) is bounded on a slightly smaller neighborhood of {a{t) : t G [0, T] }. Since 
\\{x - a{t))'^(j)j{A{t),B{t),h,a{t),ri{t),x)\\ has order ^''^'/^ it follows that ||^(-,t,^)|| has 
order /i'^^. Applying Lemma 2.8 of |[l^, we learn that the ilj{x,t,h) solves the Schrodinger 
equation up to an error whose norm is bounded by Ci h!^'^ , when U is sufficiently small. 

Note that the argument above requires the insertion of cutoffs to handle the Gaussian 
tails or some other assumption, such as 1/ G C'+2(iR^) with | D'^Vix) \ < M„ exp(rx^) for 
\m\<l + 2. 
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5 Estimates of the Expansion Coefficients 



In this section we study the behavior of Ck(t). 

The first step is to get a good estimate of the operator norm of the bounded operator 
'P|j|<n, where -P|j|<n denotes the projection onto the span of the (pj with \j\ < n. 



Lemma 5.1 In d dimensions, 

(X - a)™ P\j\<n = -P|i|<n+|m| {x - tt)™ i^j|<n, (5.1) 

and 

||(x-ar/^,,„|| < (y2AdM||)'"' {^-^^y. (5.2) 



Proof: Formula (2.22) of ^ states that 

(xi-tti) = Aip^p(A,5, ;i,a,r7)* + ^ v4ip^p(v4,5,^,a,?7)^ . 

Note that the right hand side contains 2d terms. Suppose v is any vector in the range of 
P\j\<n- Then using formulas (2.8) and (2.9) of DTOf, we easily deduce that 

\\ Ap{A, B,h, a,!])* v\\ < ^/n + 1 \\v \\, 

\\ Ap{A, B,h,a,ri) v\\ < v^ll'^^ll) 

and that both Ap{A, B, h, a, rf)* v and Ap{A, B, h, a, rj) v belong to the range of -P|j|<n+i- 
It follows immediately that 



\\{xi- tti) P\j\<n\\ < v2h d\\A\\y/n + 1, 

and that {xi - a^) P\j\<n = -Pb|<n+i - a) P\j\<n- 

The lemma follows from these two results by a simple induction. H 

f k~l\ 

The conclusion to the next lemma contains the binomial coefficients , . For 

\p-l J 

k = 1 and p = 1 we define this to be 1. 



Lemma 5.2 Suppose V satisfies the hypotheses of Theorem \3.][ . 

Fix T and choose a classical orbit a{t) for <t <T. The hypotheses guarantee that 

D, = j 1, =up f"! ^AE^XMM ] ,5.3) 

I o<|n|, o<i<r nl J 
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and 

^ " ^ (5.4) 



D2 = max J 1, sup V2 d6'^ M(^) II \ 

[ 0<t<T J 



are finite. 

We define = ( ^ ~'~ ^ I ' which is the number of multi-indices m with \m\ = 3 in d 



dimensions. 

Suppose Co(0) is a normalized vector with Cqj{0) non-zero only for \j\ < J, and suppose 
Cfcj(O) = for all j when k > 1. Let Ckjit) be the solution to U-^ )) ■ ■ ■ > M-^V ) wi-th 

these initial conditions. Then for t G [0, T], we have 

(^o,j(t) = whenever \j\ > J, (5.5) 

||co(t)|| < Di, (5.6) 
and for k > 1, 

Ck{t) = E ct\t), (5.7) 
p=i 

where 

CfciiW = whenever \j\>J + k + 2p, (5.8) 

and 



Proof: The finiteness of Di and D2 is standard. 



The conclusions ( |5.5| ) and (|5.6| ) are trivial. 

We assume t G [0, T], and let X{t) denote the formal vector whose entries Xi{t) denote the 
infinite matrix that represents h~^^'^ (xj — ai(t)) in the basis { (j)j{A(t), B(t), h, a(t),ri(t), ■) }. 
From ( [4.9|) we have 

tc.it) = K,{t)c,{t) = Y: (-P"^^)["(^)) x(t)-Co(t). 

\m\=3 

We integrate to obtain ci{t) = c^i\t). Lemma (|5.5| ), and (|5.6| ) imply two conclusions: 
cS^l.(t) = whenever \j\>J + 3, (5.10) 



and 



c™(t) II < DiDlDs 



1/2 



(5.11) 
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where the factor ( ~'~ ^ | is the number of muhi- indices m with Iml = 3 in dimensions. 
\d-l J 

This proves and Q for = 1. 

For k = 2, we have from (|4.1CI|) , 

^ ^ m^x(trco(t)+ E ra^x(trc.(t). 

|m|=4 |m|=3 

The two terms on the right hand side of this equation produce two terms, C2^^(t) and C2^^(t), 
when we integrate to obtain C2(t). Using (p.5[), ( [5 .61 ), ( |5.10| ), ( |5.11| ), and two apphcations of 
Lemma we learn that C2(t) = C2^'(t) + C2^k^), where 

cW.(t) = whenever |j|>J + 4, (5.12) 
cg(t) = whenever |j|>J + 6, (5.13) 



rf-1 —/T 



1/2 



and 



cf{t)\\ < DiDlDl y^^^y (5.15) 
(The factor of ^ ^ ~'~ ^ ^ in (|5.14|) is the number of multi-indices m with \m\ = 4.) This 



d + 3 \ ^ f d + 2 



imphes (|3), (|3), and (I]!) for A; = 2 because ^ ^ _ j < =^3- This 

combinatorial inequality follows because d > 1 implies 

f d + 2Y f d + 3Y' (rf + 2)!2 4! (rf-1)! 



d-1 J \d-l J 3!2 (d- 1)!2 (d + 3)! 

4d d + 2 d+1 



> > 1. (5.16) 

d + 3 3 2 - 1 + 3/d - ^ ^ 

From ( [4.11| ) with = 3, we have 

tC3{t) = K5{t)co{t) + K^{t)ci{t) + K^{t)c2{t) 

fc=0 |m|=5-fc 

Using d]!), ( pJOD , ( |5TTD , ([5J2| ), ([57[3| ), (|57[4l ), ( pJ[5| ), and four applications of Lemma 

Owe learn that C3(t) = C3 ' (t) + cf^(t) + C3^'(t), where 



cgl.(t) = whenever |j|>J + 5, (5.17) 
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'^f\it) = whenever 



[31 

C3j(^) = whenever \j\ > J + 9, 



> J + 7, 



4%) II < {'/_1]d,di(^'-'"^'- 



1/2 



J! 



(5.18) 
(5.19) 

(5.20) 
(5.21) 



and 



c^\t)\\ < DlDlD' 



/(j+9)!y/' i_ 

^' J! j 3!' 



(5.22) 



This imphes (pT7|), ( |5.8| ), and ( |5.9| ) for A; = 3 because of ( p.l6| ) and the similar inequality 

3 



d + A \ d + 2 

d-1 ] - I d-1 



This inequality follows because d > 1 implies 

5 



d + 2Y ( d + A \ ^ 



d+2 d+l 



d- 1 



d- 1 



l + A/d 1 + 3/rf 



> 1. 



Now suppose inductively that the lemma is true for all k < q, for some q > 2. By 
integrating ( ^.11| ) with n = g + 1, we can decompose 



<? n+l 

c,+i(t) = cWi(t) + E E d[q,n,p]{t), 

n=l p=2 



where 



4^jl(t) = -W i^g+3(-s)co(s)rfs, 


d[q,n,p]{t) = -i i^g+3-n(s) c|f~^](s) rfs, 
Jo 

for 1 < n < g and 2 < p < n + 1. We interchange the sums in ( p.23| ) to obtain 
c.+i(t) = ,(t) + E 4t,(t), 



(5.23) 

(5.24) 
(5.25) 



p=2 



where 



n=p— 1 



(5.26) 



(5.27) 
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for 2 < p < g + 1. This establishes ( |5.7| ) for k = q + 1. 

The induction hypotheses, formulas ( [4 .61) , ( |5.24|) , (|5.25|) , ( |5.27|) , and Lemma |5]l| imply 
(5.8) for = g + 1, as well as the two inequalities 



d + q + 2 
d- 1 



1/2 



J! 



(5.2^ 



and 



d + q + 2 
d 



2 - n \ +3-n / {J + q + l + 2p)\ \ 

1 j ^^^^ \{J + n + 2p-2\ ) 



1/2 



n — 1 
p-2 



{J + n + 2p-2)\ 
J! 



1/2 



n — 1 
p-2 



d+q+2-n 
d-l 



(J + g + 2p+l)!y/' 



X 



J! 



(5.29) 



From these inequalities and (|5.26| ), we obtain (|5.9| ) for A; = g + 1 as soon as we establish both 
the inequality 



d+q+2-n\ ( d+2 
d-l ] \ d-l 



< 



d + 2 
d-l 



9+1 



for < n < q, and the identity 



p — I 



E 

n=p— 1 



n — I 
p-2 



(5.30) 



(5.31) 



for q >2 and 2 < p < q + I. 

We set r = g — n and note that ( |5.30| ) is equivalent to 

-r-l 

< 1, 



d + r + 2 
d- I 



d + 2 
d- I 



(5.32) 



for < r < g. However, 

/ d+r+2' 
[ d-l ^ 

(rf + 2 + r) 



-r~l 



d + 2 
d-l) 

/ {d-l)\ 3! 

{d-iy. (r + 3)! 1^ (rf + 2)! 



r+l 



r + 4 /1\ ''+^ 
,4, 



r + 5 /2\''+i 
.5, 



r + d + 2 fd-V 
d-l [d+2 



r+l 
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Inequality ( p.32D follows if each of the factors in the square brackets is bounded by 1. Thus, 



we need only prove 

r + m + 3 ^ f m + 3 ^'^^ 



m \ m 

for 1 < m, which can be verified by using the binomial expansion: 

r + m + 3 ^ r + 3^^ , ^,3 / 3^+^ / m + 3 

= 1 + < 1 + r + 1 — + ...= 1 + — = 

mm m \ mj \ m J 

This proves (|5.32| ) and hence ( p.30| ). 

The identity (|5.31|) is trivial for g = 2 and p = 2, 3. Assume inductively that it is true for 
all 2 < g < m and 2 < p < g + 1, where m > 2. The identity ( p.31| ) is trivial for q = m + 1 
and p — 1 = m + 1, since 

/m + l\__|^_/m\ 
Im + ly \ m I ' 

Then for m + 1 > p — 1, we have 

f m + 1 \ (m + 1)! 



p-1 J {p-l)\ {m-p + 2)\ 



{p — 1) (ml) ^ {m — p + 2) (ml) 



{p-l)\ {m-p + 2)\ {p-l)\ {m-p + 2)\ 

I m \ I m \ 
[p-2 ) + [p-1 ) 



m+l 

E 

n=p— 1 



n — 1 
p-2 



This proves (|5.31| ) and completes the proof of the lemma. 



Corollary 5.3 Assume the hypotheses of Lemma \5.2^ . Then in addition to ( |5.5y and ( |5. 6| j, 

we have the following for k > 1: 



Ckj{t) = 0, whenever |j|>J + 3A;, 



(5.33) 



and 



Ck{t)\\ < 



(J + 3fc)! DiDl 



k\ 



1 + DiDU 



(5.34) 



Proof: Since p < k, ( pTSD implies ( |5.33| ). 
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To prove ( |5.34|) , we note that j is increasing in p. Thus, (|5.9| ) and p < k 



imply 



Summing over p, we obtain 



This imphes ( ^.34|) . H 

6 Optimal Truncation Estimates 



In this section we show that the error given by ( [4.12| ) and (|4.13| ) is exponentially small if we 
choose Z = [(7/^] for an appropriate value of g. 

The philosophy will be separately to estimate the error near the classical orbit and far 
from the orbit. To do so, we let h be any positive number and define Xi(2^5^) to be the 
characteristic function of { x : |a; — a(t)| < 6}. We set x^i,^-,'^) = 1 — 



Lemma 6.1 Assume V satisfies the hypotheses of Theorem p. i| . Define Xii^jt) as above 
and ^i{x,t,h) by ( U-lSj ). For fixed T > and b > 0, there exists Gi > 0, such that for 
each g G (0, Gi), there exist Ci and 71 > 0, such that if I is chosen to depend on h as 
1(h) = Ig/h], then 



rT 



/ \\ Xii; t)^nn)i;t,h)\\dt < C,exp{-^,/h}. {Q.l] 







Proof: It is sufficient to prove the existence of ai and jSi, such that 



T 



/ \\ Xii; t)U;t,h)\\dt < a^l3[l"^h"\ (6.2) 







If this can be established, we choose Gi = /Jf^. Then < g < Gi and / = [(7/^1] imply 
f3f g = e~'^, with u > 0. Since ai{f3flhy/'^ = ai e~'^^^^'^^\ this implies the lemma with 
Gi = ai and 71 = ujg/2. 

To prove ( |6.2| ), we note first that our hypotheses imply the finiteness of 

=up ^Nl(g-''')WI. (6.3) 

|n|>0, 0<t<T, \x-a(t)\<b n\ 
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We use this, (|]T3|), and to see that 

\\Xi{x,t)^i{x,t,h) II 



< 



l-l 



k=0 \m\=l+2-k 



ml 

X E Ckj{t)<Pj{A{t),B{t),h,a{t),r]{t),x) 
\j\<Jii) 



k=0 



2+k 



X E 

\m\=l+2-k 



{x-a{t)r E CkAt)MA{t),B{t),h,a{t),v{t),x) 

\3\<J{1) 



l-l 



< D,h"'+' E E \\x{trc,{t)\\, 

k=0 \m\=l+2-k 

where X{t) is the infinite matrix that represents h~^^'^{x — a{t)) in the (pj basis. 
Thus, 



h-' r \\xi{;t)U;t,h)\\dt < D,h"'Y. r 6-'-'+' E \\Xitrckit)\\dt. 

fc=0 \m\=l+2-k 

(6.4) 



We apply Lemmas |5.1| and |5.2| to estimate each integral on the right hand side of (|6 
For /c = 0, we obtain 



d + l + l\ / (J + / + 2)! 



1/2 



J! 



< 



J\ 



(/-I)! 



Di Dl T. (6.5) 



In the last step, we have used -D2 > 1, (|5.30|) , and (J + / + 2)! < 777— ^—^r^, which is true 

{(i-iy.y 

for / > 1. 

For /c > 1, we write the integral on the right hand side of ( |6.4[ ) as a sum of k terms by 
employing (|5.71 ). By (|5.8| ), ( ^.9|) , and Lemma ^TT], the p"^ integrand satisfies 



^-l-2+k ^ X{t)'^ct\t) 
\m\=l+2-k 



< D 



l+2-k 



d + l + l~k\ f {J + l + 2p + 2) 



1/2 



d- 1 



1^ {J + k + 2p)\ 
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< 



k-1 
p — 1 

k-1 
p — 1 



Dl D 



J\ 



d + l + l-k\ / (J + / + 2p + 2)! 



J! 



p\ 

p\ 



(J + / + 2p + 2)! 
J! 



1/2 



(6.6) 



In the last step, we have again used ( |5.30|) . 

We now mimic the proof of Corollary p.3| and then integrate to obtain the following 
estimate of the k^^ term in (16? 



(J + 3/)! 
J! 



(/-I)! (A; + l)DiD^ 



1 + DiD^T 



k+l 



- 1 



We define 



1 + ^1^2'^, 



bound 



D 



{k + l)D,Dl 



by D^, and [(1 + DiD^T) 



fc+i 



(6.7) 



by -Ds^^. We then sum over k in 



to obtain the estimate 

■T 



-1 







■,t)6(-,t,^)Nt < /^4^'/' 



(J + 3/)! 
J! 



(/-I)! -1 
by 



(6.8) 
— , which 



In this expression we bound (J + 3/)! by (J + 3/)"^+^' and n, 

(/-I)! p'-i(/-l) 

holds for some constant p. After some algebra, this leads to the estimate ( |6.2|) . I 

To prove the analogous lemma with xi replaced by X2, we use spherical coordinates for 
In spherical coordinates when d > 2, the operator — Aj^ + has the form 



d-ld 



£2 



r dr 



7^ + ^ + r^ 



Here C!^ is the Laplace-Beltrami operator on S*"^ ^. For d > 3, it has eigenvalues 



qiq + d-2), 



(6.9) 



with multiplicities 



(rf-1) 
1 

Jd^. 



^(g + l)(g + 2) ■■■ (g + d-3) | (g + d - 2)(g + d - 1) - (g - l)g | 



+ 1) (g + 2) ■ ■ ■ (g + d - 3) (rf - l)(2g + - 2) 



(6.10) 
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where q = 0, 1, . . .. We denote a corresponding orthonormal basis of eigenf unctions by 



Yq,nii^) for 1 < m < niq. 



even functions and multiplication by A = 1 on odd functions. The operator — |^ + on 
IR just becomes the direct sum of two copies of — |^ + on (0, oo) with Neumann and 



When d = 1, the analog of the Laplace-Beltrami operator is multiplication by A = on 

Dirichlet boundary conditions at r = 0. 

When d = 2, ( |6.1UD should be replaced with niQ = 1 and nig = 2 for g > 0, but the 
inequality rriq < Cde"''''^ still holds. 

The eigenvalues of — A„ + y"^ are E = An + 2q + d with normalized eigenf unctions 



Here 



\r{q + n + 



r« Ll^'^~\r') e-^'/' Y,,^. (6.11) 



2' 



„m 



Li(.) = t i-ir ^ (6.12) 

m={) \ / 

denotes the Laguerre polynomial that satisfies the differential equation 

xu"{x) + {[3 — X + 1) u {x) + nu{x) = 0, (6.13) 
and the normalization condition 

r L^r^ix) Liix) x^ e-^ dx = I r(/n+i) |! ^ , (6.14) 

for /? > -1. 

The following lemma implies an estimate for \il)q^n,m{f',i^)\ when r is in the region which 
is classically forbidden because of energy considerations. 



Lemma 6.2 For /5 = g + | — 1 with q = 0, 1, . . ., the Laguerre polynomial L^{x) in (\6. 1 ]\) 



satisfies \L^{x)\ < — whenever x > An + 2f3 + 2 = 4n + 2q + d. 
n] 



Proof: We mimic the proof of Lemma 3.1 of |[TT|. The first step is to show that g{r) = 
L^{r'^) cannot vanish in the classically forbidden region > An + 2q + d. This 

function vanishes at infinity and is a non-trivial solution to an equation of the form 

-g"{r) + w{r)g{r) = 0, 

where w{r) > for > 4n + 2g + d. From this differential equation we conclude that g 
and g" have the same sign in this region. By standard uniqueness theorems, g and g' cannot 
both vanish at the same point. To obtain a contradition, suppose g has a zero at some point 
ri with rf > 4n + 2g + d. Since g vanishes at infinity, the mean value theorem guarantees 
that g'{r2) = for some r2 > ri. Without loss of generality, we may assume g{r2) > 0. This 
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forces g"{r2) > 0, so g' is locally increasing. It follows that g' is increasing for all r > r2. 
Thus, g could not go to zero at infinity. This contradiction shows that g could not have had 
a zero in the region. 

We now proceed by induction on n. Since Lq{x) = 1, the lemma is true for n = 0. We 
now assume n > 1 and that the lemma has been established for L^_i{x). 
Our non- vanishing result and (|6.12|) imply 

= 1 - B,4^), 
(-1)" ^ ' 

where Bp^nix) = 0{l/x) for large x, and B[^^n{x) > —1, for x > 4n + 2/5 + 2. 
Using recurrence relation 8.971.3 of we have 

d ^ , , xL'l'ix) — nLf^ix) n\ 



dx ^ ' (-1)" 

^n+l (-1)"' 

By our induction hypothesis, Ln-i{x) has sign (—1)""^ for x > An + 2(3 — 2, which includes 
the region of interest. 

Thus, Bjs^nix) is increasing. Since it goes to zero at infinity, it cannot be positive. This 
implies the lemma. H 



Lemma 6.3 Assume V satisfies the hypotheses of Theorem 3_A. Define X2{.x,t) as above 



and ^i{x,t,h) by For fixed T > and b > 0, there exists G2 > 0, such that for 



each g G (0, G2), there exist C2 and 72 > 0, such that if I is chosen to depend on h as 
l{h) = Ig/h], and h is sufficiently small, then 

r \\X2i;t)^i^n)i;t,h)\\dt < C2 exp { -72 } . (6.15) 



Proof: We begin by using the analyticity of V to control Taylor series errors. We define 

Csix) = {ze(D'^ : Zj = Xj + 5e*^^ 9j E [0, 2n), j = 1, 2, ■ ■ ■ , d}. 
If z E Cs{C{^^ ct)), then, for all j = 1,2, ■ ■ ■ ,d, 

\zj\ < 6+\(j{x,a)\ < 6 + \aj\ + \xj — aj\. 
Using this and applying < 2(6^ + c^) several times, we see that z G Cs{C{^: C')) implies 

I V{z) I < M exp(2r(x - a)^) exp{4T{6^d + a^)). 

Hence, writing — -D'^V{C{ X ^ Qjj J cLS cL d-dimensional Cauchy integral, we obtain the bound 
m! 
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where ({x, a) is any value between x and a. Thus, for < t < T, there exists a constant 
Ml, such that 



l|(W)(ax,a(t)))| < ^ exp(2r(x-a(t)f). 
We use this, j ^J^ ), ( ^1^ , and to see that 

II X2ix, t)^i{x,t,n) II 



(6.16) 



< 



i-i 

E 

fc=0 



|p|=i+2-A: 



\p\=l+2-k 



k=0 



pi 



X E Ck,j{t)(Pj{A{t),B{t),h,a{t),r]{t),x) 

\j\<J+3k 

X2{x,t) exp(2T(x - a{t)f) (x - a{t)Y 
X 5] Cfe,,(t)0,(v4(t),5(t),;i,a(t),r/(t),a;) 

|i|<J+3fc 



l-l 



|j|<J+3fc |p|=«+2-fc 



fc=0 



X (/),(A(t),S(t),;i,a(t),r/(t),x) 



(6.17) 



(6.18) 



Note that (|5.33|) has been used to replace |j| < J{1) with |j| < J + 3/c. 
The norm in the final expression of (|6.17|) equals 

X{z: \z\>b}{x) exp(2rx^) x^ 0i(A(t), B{t), h, 0, 0, x) 

where |p| = I + 2 — k. 

We assume that h is sufficiently small that 4r/i|A(t)p < 2/3. Then the square of the 
quantity ( |6.18| ) equals 

2-lil vr-'^/2 I det A(t)|-i ^^^^/^ 

X / a;2V™'|H,(A; |A(t)ri ;i-i/2^)|2 e-I^WI"'^'/^ 



|x|>fe 



2\\P\ 



< 



< 



2lil(j!)7r'^/2 Jn^/^\\A(t)\y\>b 

in\\Ait)rf 



< e 



2lil(j!)7r'^/2 Jn^/^\\A{t)\y\>b 

(^pwin'"' 



2lil(j!)7r'^/2 7ni/2||A(t)|j;|>b 
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|ypH|7^,(A; y)\'e~y'yUy 

\y\'^P^\nM;y)\'e~'"'^'dy. (6.19) 



By formula (3.7) of [TlJ, nj{y) 



1 



Hj{A; y) e ^^^^ is a normalized eigen- 



function of — + with eigenvalue 2\j\ + d. Thus, in spherical coordinates, it can be 
written as 



{q,n,m: 2n+q=\j\} 



1. 



(6.20) 



where ^ Mi,g,n,mP 

{i},n,m:2n+(j=|j|} ^ 

We ultimately choose I = {g/h], with < g < G2. Since J(/) = J + 31 — 3, there exists 
C3, such that ^ < 1 implies J(/) < C^/h. By choosing 6*2 sufficiently small, we also have 
J{1) < {\\A{t)\\'%'^ - l)/{2h) for < t < T and small h. Thus, the relevant values of j in 
( Fl9| ) satisfy J2\j\ + d < \\A{t)\\-^bh-^^\ 



Lemma |6^ shows that 



So, we see that ( |6.19| ) is bounded by 



< — — whenever r > 4n + 2q + d = 2\j \ + d. 



2bl 



||A||r1/2 



{q,n,m: 2n+q=\j\} 



,rV2 r'i-^drdu. 



We interchange the sum and integrals and apply the Schwartz inequality to the sum. This 
shows that ( |6.19| ) is bounded by 

2n\ 



^ T o -\- - -\-n) 

{q,n,m:2n+q=\j\} ^ ^ 2 ^ 



r 



d-l+2|p|+2q 



-bV(6||A(t)||2n) ^bi p(^)||2bl y 



Lf^-Vr^) 



'^''^dr 



2n\ 



{q,n,m: 2n+q=\j\} 



X 



r 



d-l+2\p\+2q 



A||fil/2 



(6.21) 



By reducing the value of G2 if necessary, we can ensure that the hypotheses of Lemma |6.2| 
are satisfied in the integration region in the right hand side of ( |6.21j ). So, Lemma |6.2| shows 
that the integral satisfies 



^d-l+2\p\+2q 



\\A\\h 



-'"'^dr 



< 



[nl 



^4n+d-l+2\p\+2q g-r-V2 
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(n!)^ Jo 

22n+f+g+|p|-l ^ 

= 7^, T{2n + - + q+\p\). 



So, (|6.21|) is bounded by 



{g,n,m:2n+<7=|i|} ^' + ^ + 

We use ( |6.10| ) to estimate the sum over m < and bound this by 



^ -62/(6||A(t)|pn)^H||^..N||2b| a,, 2l-^l+2+H r(|j| + f + IpI) 

{n...t^|-2n} ^! r(g + f + n) 

C,e-^'/(6ll^WII'^) /ilPl IPl 2l^-|+Mfl e-'^l^'l r(|j| + ^ + 



g-2«d" 



Since e"^""^" < 1, this is bounded by 

2 „<|,-|/2 ^! r(g + f + n) 

(6.22) 

For n < \j\/2 and fixed, there exists a constant C", such that r(|j| —n+—) > C" (|j| —n)\. 
So, the sum over n in ( |6.22 ) is bounded by 



Thus, ( |6.21| ) is bounded by 



^„„^_fe2/(6p(t)l|2^) ^bl p(t)f 1^*1 21^*1 e^^l^'l ^^^^^ ^-^^^ (6.23) 



This quantity bounds ( |6.19| ), which, in turn, bounds the square of ( |6.18| ). Terms of the form 



( |6.18| ) occur in ( |6.171 ). Putting this all together, we see that (|6.17|) is bounded by 



^^^/„/l/2^_fe2/(i2j|A(t)|l2a)^-«-2 ^ ^fc/2^fc |Cfc,i(^)l 1.7 1/2 



\\A{t)t\ 2H/2 Jr(|j| + ^ + bl) . (6.24) 

|p|=i+2-fc ^ 
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The number of terms that occur in the final sum of this expression is| 



the terms in that sum are increasing. Thus, ( |6.24| ) is bounded by 



i-i 



^==0 lil<J+3fc v|j|! 



^,i^^.„„i/2^-iP,ii2,Ai,)rk,g-,-2 (2n\\A(t)fy*' 'y, {2\\A(t)\f 



k=0 



d-1 

Applying the Schwartz inequality to the sum over j, we see that this expression is bounded 
by 

(\ 1/2 

The number of terms that occur in the final sum of this expression is ~^ '^^ ~^ ^ ^ , and 
the terms in that sum are increasing. Thus, the expression is bounded by 



j^^^..i/2^-by(i2\m)fn)^-i-2 ^2h\\Ait)fY Yl 



fc=0 



/ l-k + d + l\ 

M d-1 11^^ 



d J {J + 3k)\ ^ 2 
We now apply the estimate of ||cA;(t)|| from Corollary |5.3| to bound this by 



,/3,(j+3fc)/2 + D,Dl\t\r (T{J + 2k + 1 + 1 + 2 



'"^ -k/2 f l-k + d+l\ f J + 3k + dY^^ 

[^\\^[^)\r) 

k=0 \ 



d — 1 I \ d I 



k\ \ J\ 



ill cPdJ/2 



J! 
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k=0 



d — 1 \ d 



X e^^^^/^ \\ \T{J + 2k + - + l + 2). 

rv. V 2 

(6.25) 

We now employ the following inequalities that hold for some numbers and Dq: 
1^ l~k + d + l ^ < 1^ l + d + 1 ^ 

( J + 3k + d\ (j + 3l-3 + d\ 
[ d ) - [ d j' 

D',Dl{l + D,Dl\t\f < Di-\ 



l_^nj + 2k + ^ + l + 2) < Jf^fiJ + ^l + ^)^ and 

k=0 

We then see that ( |6.25D is bounded by 



2\ 2 



X 



(6.26) 

We bound this expression by using the two inequalities 



/ l+d+1 
[ d-1 

J + 3l-3 + d 
d 



< {l + d+lf-' and 



< {J + 3l-3 + dy 



Note that the right hand sides of these inequalities grow polynomially with /. 

Since d is fixed, we conclude that (|6.26|) , and hence, (|6.17|) are bounded by a constant 
times 



^ (/-I)! V ' 2' 
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for some positive 7, and 7'. With J fixed, we apply Stirling's formula to the factorial and F 
function to bound this by another constant times 



^-H' /{12\\ A{t)W'h) ^1/2 (30^'^^ 



By choosing I = {g/hj for some sufficiently small g > 0, this is bounded by a constant 
times e~^^^^. 

This implies the lemma. H 



Theorem |3.1| follows immediately from Lemmas (|6.1|) and (p. 3D with G = min { Gi, G2 } 



7 Localization Estimates for the Wave Packets 

In this section we show that our wave packets are localized near the classical path. Given any 
e > 0, we can choose the truncation parameter > 0, such that our exponentially accurate 
wave packet is concentrated within {x : \x — a{t)\ < 6 } up to an exponentially small error. 



Proof of Theorem |3.2| : Let x{^i'^) be the characteristic function of the set { x : |x — 



a(t)| > 6 }, and let ^/'(x, t, fi) be the result of our construction with the series truncated with 

m = ig/ni 

We must prove 

\\x{;t)tlj{-,t,h)\\ < exp{-r/;i}, (7.1) 

for some F > when g > is sufficiently small. 
The left hand side of ( [7. 1|) is bounded by 



E E Mt)\ \\xi;t)UMt),m,^^{t),vit),-)l (7.2) 

^==0 \j\<J+3k 



The norm in this sum has the form (|6.18|) , with n = and r = 0. Mimicking the estimation 



of ( |6.18| ), we obtain the estimate that corresponds to ( |6.24| ). We conclude that if (7 > is 
sufficiently small, then 

\\xi;t)UAit),Bit),h,ait),vit),.) \\ < e-^Vd^il^wip.) 

whenever |j| < J + 3l{h) — 3, for some Pd- 

We use this and the Schwarz inequality to obtain, for some constants Co, Ci, C2 and C3 

E E KM \\x{;t)4>Mit),B{t),h,a{t),7^{t),-)\\ 

^=0 |j|<J+3fc 

l-l ( \ 

< E II Cfe(t) II e"^'/(i2||AWipn) ^ ^\m, 

/c=0 \ |j|<J+3fc 
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fc=0 

oo 



A;=0 

provided g is small enough that hC2k < €29 < 1 is satisfied. H 

8 Ehrenfest Time Scale 

In this section we consider the accuracy of our construction when we allow T to grow as 
^ — 0. Since the results stated in Theorem p.4| and the method of proof are basically equiv- 



alent to those of |[TT|, we will be rather sketchy. 



Proof of Theorem |3.3| : The first point to notice is that the since potential is bounded 
from below, energy conservation implies that a{t) grows at most linearly with time. The 
exponential bound on the potential then implies the existence of Z)i > and f > such 



that the quantity (|5.3|) is bounded by 

Di{T) < 1^1 e^^^. 

Similarly, the existence of the Lyapunov exponent A implies the existence of D2 such that 
the quantity ( ^.4| ) satisfies 

D2{T) < Die^^. 

It then remains for us to keep track of the time dependence in the proof of Theorem 



In particular, the quantities ( |6.3D and (|6.7|) fulfill the following estimates, modulo a possible 
increase of v. 

D,{T) < D,{T), 
D^{T) < D5e(^^+^^)^. 

Using these bounds, we get the existence of constants C and independent of time, such 
that ( |6.2|) can be replaced by 



Thus, if we choose I = g(T)/h, then (|8.1| ) is bounded by 
SO that we need 

g{T) e(6^+2-)'r _ and g{T)/n^ 00. 
These demands are satisfied by the choices 

g{T) = e-^T and T = T'ln(l/;;,), (8.2 



1/2 



.11 
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provided 

6\ + 2vT < K < l/r. (8.3) 

Note that the prefactor in will be of order h^'^^ , for some finite Ui. It will thus play no 
role since it follows from these considerations that there exists 71 > 0, such that 



T 



/ \\xi{-,tM-Xmdt = 0(^-^^e-^^/^ 



1-kT' , 



By a similar argument, we obtain the estimate corresponding to ( |6.15| ) with other time- 
independent constants D' and C: 

h~'J^ \\X2{;tM;t,h)\\dt < h-'e-''/^''^'"'^^D' (C'/^e(«^+2..)Ty/2_ ^^^^^ 
Inserting our choices (|8.2|) and constraints (|8.3|) in (^.41), it is elementary to see that there 



exists positive 1^2 and 72, such that 

Jo 

which proves the Theorem. I 



Proof of Theorem 3.4: Considerations similar to those in the second part of the proof 



of Theorem p.3| show that there exists constants Co, Ci independent of T such that 
\\X2{;t)ij{-,t, h)\\ 

l-l / \ 

< E^'^'Mm E \\X2{;t)UA{t),B{t),n,a{t),r^{t),.)r 

k=0 \ |i|< J+3fe 

fc=0 k=0 

where, by virtue of ( ^.2|) and ( ^.3|) , we can take Ci < 1. So, the Theorem holds with 
exponential decay of order q-^^ / C^"^^^^ ^^"^ | 

9 Scattering Theory 

In this section we show our approximations are valid up to exponentially small corrections 
in a scattering framework, provided the potential satisfies hypothesis D. 



Proof of Theorem |3.5| : First note that equations ( |3.2| ) together with 



for any j E N"* imply that as t ^ ±00, 

^itHo{n)/h ^iS{t)/n (j,^^A{t), B{t), h, a{t),r]{t), x) e'^^'^ 5±, h, a±, r^^, x) 
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with S- = 0, for any j G N"^. Moreover, using ( |3.2| ) and the property 

mm{\v\,l){t) < {tv) < max(|t;|, 1) (t), 

for any v G IR!^ and any t G IR, with (t) = y/l + 1^, we get the existence of Cq > and 
Ci > depending on the asymptotic data {a±,ri±), such that 



D"1/(a(t)) 



< 



~ ~|m| 
(t)/3+l"^l 



(9.i: 



for large times. This estimate together with ( p.3|) and Lemma yields the following 
estimate on the operator Kk{t) P\j\<n defined in ([4.6[): 



\Kk{t) R,|<„| 



< 



\m\=k 

d-l + k 
d-1 



'{n + k)\ Co 4 
n\ Iff 



(9.2) 



where £2 depends on the asymptotic data {a±,ri±, A±, B±) and the binomial coefficient gives 
the number of multi-indices of order k. At the possible cost of an increase in the constants, 
we may assume this estimate is valid for all t G IR. 

This estimate shows in particular that Kk{t) P\j\<n is integrable in time. From this, it 
is easy to check inductively that the solutions c„(t) to the equations ( 4.11| ) have limits as 



cxo. 



The asymptotic values of the coefficients c„(t) at infinity allow us to define the asymp- 
totic states ^±{A±, B±,h,a±,ri±,x) by ( p.6[ ) with initial conditions at —00 characterized by 
arbitrary normalized coefficients that satisfy 



Co J (-00) 



0, for \j\ > J, and 

0, for n = 1, 2, ■ ■ ■ , and all j G N"^. 



(9.3) 



Thus, our approximate solution 

'ilj{x,t,h) = e*^W/^ c,{t,h)(P,{A{t),B{t),h,a{t),r]{t),x) 

\j\<J+3g/h-3 

has the asymptotic property as t ^ ±00, 

e^*^o(^)/^V^(x,t,;i) ^ $±(A±,i?±,^,a±,r/±,x). 
We prove below that 

(x,t,n) - jim e'(*-^)^(^)/^e'^^«('')/^$_(A_,5_,;i,a_,r/_,x)|| 

[x,t,n) - e^*^(^)/^^]+(/^)$_(A_,5_,;^,a_,r7_,x)|| 
0(e-^/^), 



(9.4) 
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uniformly for t G IR. Thus, making use of ( |9.4| ), we have 

hm ||e'*^»(^)/^7/'(x,t,^) - e**^«(^)/^e'*^(^)/^l]+(/i)<l>_(A_,5_,;i,a_,r7_,a;)|| 

t — ^'^--CxD 

= ||$+(v4+,5+,^,a+,?7+,x) - ^l~{h)* n^{h) ^_{A_, B^,h,a_,r]^,x)\\ 
= \\<^+{A+,B+,h,a+,r]+,x) - S{h) <i>^{A^, B^,h,a^,r]^,x)\\ 
= 0(e-^/^). 

Hence, we need only show that the estimate on C,i{x,t,h) corresponding to our approxi- 
mation yields an exponentially small correction term after choosing I = g/h for sufficiently 
small uniformly for t G IR. 

We mimic Section 5 to get estimates on the coefficients 

Ck{t) = jz c¥(t) (9.5) 
p=i 

starting with 

co,j(t) = coj(-oo), coj(-oo) = if |j| > J, and ||co(-oo)|| = 1. 

We note that the number of components of the vectors c^'(t) is the same as in ( |5.8|) and 
that the combinatorics associated with the n and p dependence of the estimates is identical 
to that performed in Section 5. 

Hence, with D^^ = ^ ^ ^ ' ^^^^ order we have 

llciWII = \\c^\t)\\ < Ds^f^^~Co~4 f_J^)-'ds. 
At second order, we obtain C2(t) = C2^'(t) + (^2\^)i where 

llcf'WII < Dl cocl f Js)-P ds and 

llc^WII < Dl ^^-^-^ clcl ds, {s,r' f_l ds, {s,)-^. 
At third order, we obtain C3(t) = (t) + c§\t) + c's^t), where 



-00 



\\c?\m < 2D? 



3 



f ds,{s,)-^ r ds2{s2)-^ and 

J '. J —CO J —CO 



II '-3 



< ^3' f ds,{s,r^ r ds,{s,)-^ r dssis,)-^. 



> —CO J —CO 
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Using the identity 

-oo J — co J —CO Tl\ \J— oo 

we get estimates identical to (|5Tl| ), (|5l^ ), (|5l5|) , (|5:20|) , ( ^ ), ( ^ ). 

It is easy to check that the induction can be carried out exactly as in Section 5 to give 
an analog of Corollary 5.3 that states 



Lemma 9.1 Assume the decay hypothesis D. The expansion coefficients (\9.3^ ) satisfying 
( \9.Si ) obey the following estimates: 

Ckj(t) = 0, whenever |j|>J + 3fc, 



M 



for p < k, and 



Ck 



Our next task is to estimate the norm of ^i{x, t). We again consider separately the errors 
near the classical orbit and those far from the orbit. Let b{t) be a real valued function that 
satisfies 



(9,6) 



for all t E IR. We define Xii^^t) to be the characteristic function of {x : \x — a(t)| < b{t)} 
and X2{x,t) = 1 — Xii^^t)- Then, for some constants £3 and £4 and any t G IR, we have 



Xiix,t)) 



< 



< 



\m\ 



(C„(a;,a(t))))/^+IH 

~ ~\m\ 
C3 Q 



Xiix,t) 



(t)/3+IH ' 

since for large times on the support of Xi, 

\(m{x,a{t)))\ > |a(t)|/4. 
Therefore, for some constants C5 and Cg, 



Xl\.X,t) — j X[t) f\j\<n 



7. 



< 



(n + |m|)! C5 Cg 
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for any t G IR. We now mimic the manipulations performed in Section 6 to get 

\\Xi{x,t)^i{x,t,h)\\ (9.7) 

E E Xlix, t) (D"'^)((rn(x,a{t)) ^1^1/2 _ ^^^^^^ ^_|^|/2 



< 



k=0 \m\=l+2-k 



ml 



X iij|<j+3fc E Ck,jit)(l)j{A{t),B{t),h,a{t),ri{t),x) 

\j\<J+3k 



< E 

A;=0 



X 



rf-l + / + 2- A; 
d- 1 



max 

{m: \m\=l+2-k} 



m! 



fc=0 



d-l + l + 2-k\ l{J + 2k + l + 2)\ C5g[i+^-*^ 
rf- 1 



J! 



X 



~fe 7-)fc 

k\ 



(l + C0C2 I ds (s) Co cl I ds (s) 

\ J— 00 / J —00 



Making use of ( |5.3U| ), the definition of D^, and introducing another constant 



C7 



1 + C0C2I) C2/C6 



where 



{sy'ds, 



( pr?!) is bounded by 

C5 Co (cg €2)^ / 



A:=0 



(J + 2A; + / + 2)! 



A;! 



(57 - 1) VJ! {t)P 

C5 Cq (c6 ^2)^ I 



(/-I)! 



(9. 



This estimate is integrable for t & IR and yields a bound that implies exponential decay in 
h by the optimal truncation technique. 

We now come to the estimate of \\x2{x,t)C,i{x,t)\\, which is a little bit more elaborate. 
The difficulty stems from a lack of sufficient information on the position of (m{x,a{t)). So, 
instead of the usual Taylor series error formula, we use the definition 

W^l{x)=V{x)- E ^^^^^^^{x-a{t)y 



|m|<q 
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We ultimately use q = I + 1 — k, where k = 0, 1, ■ ■ ■ , Z — 1. Our proof requires the space 
dimension to satisfy d > 3 in order to obtain integrability in t. 
Consider the following integral 

= I^^X2{x,tf\<P,{Aii)^B{t),^(^it),vit),x)\'V{xYdx 
\^j{A{t), B{t), h, 0, 0, 2)p V{z + a{t)f dz. 

\z\>b{t) 

We use the formula 



\(t>j{A{t),B{t)X^,Q,z)\ 



^]\2\o\T:'^/^\detA(t)\Ti'^'^ 



the asymptotic behavior ( |3.3| ), and the following estimate, which is valid on the support of 



mw - mw - i6P(t)p 

to obtain the bound 

J\A>b{t) j\2\i\T:^l^\deiA{t)\K^'^ ^ k jj , 

for some finite, positive b. Note that this estimate has a uniform exponentially decreasing 
pref actor. 

As in Section |, we use spherical coordinates and the decomposition ( |6.2q ) 

{q,n,m:2n+q=\j\} 



where ^ Mi,g,n,mP = la.nd.Vtj{y) = J ^.^ ]j^d/2 ^) ^ This leads 

{g,n,m:2n+q=\j\} ' J'^ 

to the estimate 

{,,n,„^:tr+,=b1} |detA(t)|;i'^/2 

where the spherical coordinates {rz,ujz) describe the vector h^^"^ \A{t)\~^ z. 

We choose p > 2, such that d/(3 < p < d, and define s > 2 by 1/s + 1/p = 1/2. 
Applying Holder's inequality, we get the bound 

- |detA(t)|;.'^/%,,„,„^tr+,=bi} 'VN>Kt) m,n,^y^^ j 

We need to bound the integral in this expression. We change variables to y = \A{t)\^^z/h^^'^ 
and use the estimate | |74(t)||/ 1 < b{t)/\/h, which is valid when \y\ < b{t) / {\\A{t)\\\/h) . This 
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yields 



\z\>b(t) 



< / ,,,, e+y''/'\^,,^,Ur,u;)ndetA{t)\h''/'dy, (9.9) 

J\y\> ^i) 

where the spherical coordinates (r, u;) now describe the vector y. Note that we have used 
det(|yl|) = I det(74)|, which follows from A = Ua \ A\, where Ua is unitary. 

Since 6(t)/||74(t)|| has a strictly positive infimum 6, and we ultimately choose / ~ g/Vh, 
with g arbitrarily small, we can assume the integration in ( |9.9|) is within the classically 
forbidden region where Lemma 4.2 applies, for all indices {q,n,m : 2n + q = of interest. 

Hence, manipulations similar to those performed in Section |^, show that (|9.9| ) is bounded 
above by 



Js'i-^ Jb/VK 



r(g + n + d/2)^/2n!^/2 Js-^-^ ' ' Jb/Vn 

|detA(t)|;.'^/^(2/.)'^+-l^-|2-/2 f isr(^±M\/o 

This implies the estimate 

e-^/^ \\V\\l {2/sfWs+\j\) r ((d + s|j|)/2)'/^ 



22/^-1 |detA(t)|i-2/« n'^/^-'^/' 

E 



X 



Jgd-i du \ Yq^m\ 



2/s 



y{q,n,m:2n+q=\j\} 



r(g + n + d/2)^/2n!^/2 



We bound the integral in this expression by using the following crude lemma. Its proof is at 
the end of this section. 

Lemma 9.2 For some constants Mq and Mi, we have 

in.mMI < Mo Ml (9.10) 



We use this and the inequalities q < \j\ and ( ^'^^ ) < 2^^^ to estimate 



n 



2/s 



Jsd-1 duj \ Yq 



{q,n,m:2n+q=\j\} 



r(g + n + rf/2)^/2 721^/2 



- l^^o r(|j| -n + rf/2)^/2n!^/2 
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< 



< 



2/s 



C"\j\\ 



E 

n=0 



n 



Hence, for some constants A^o and A^i, that depend on d and s only, 



\detA{t)\yph'^/^' 



By our choice of p. 



1 



|detyl(t)|Vp 
is integrable. 

For k < I ~ g/h, with sufficiently small g, this last estimate allows us to bound the cor- 
responding term in t)x2{x, t)\\ as follows (where the Ni, i = 0, 1, 2, 3 . . . are constants): 



k=0 
l-l 



J2 Ck,j{t) X2{x, t) V{x) MMt): B{t), n, a{t),rj{t),x) 

\j\<J+3k 

< Wckim [ E \\V{x)<l^,{A{t),B{t),h,a{t)Mt),x)f 

. \j\<J+3k 



k=0 
l-l 



< E 

A;=0 



^k/2 jJ + ^ky. Ni 

-h/2h 



pNo 



1/2 



< 



< 



J! A;! 



det A{t)\yPh'^/P 



E ^i'^'' 

y\j\<J+3k 



(9.11) 



j2 k^i^n^/^N^ 



(9.12) 



deiA{t)\yPh'^'P' 
It remains for us to control integrals of the form 

F\p) = X^(x, t)\ct>,{A{t), B{t), h, a{t),r]{t),x)\'\DW{amx - aimyiplfclx 



< 



< 



7.2 ~m 



2/3 



(x - a(i))2|^l |0j(A(t), S(i), ^i, a(i), ?7(i), x)|2 
l/2H|7^,-(A;y)|^e^'/V(2l^-|j!7r'^/^)dy, 



'l?/l>5/v^ 

where we used the same type of estimates as above. We bound the last integral in this 
expression by using spherical coordinates and noting that the integration region lies within 
the classically forbidden region, if g is sufficiently small. The integral is thus bounded by 

o2n+g+|p|+d/2-l 
{q,n,m:2n+q=\3\} ' ' 
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So, for some other constants we have 

p. .2 . /2ij^Vj^'^l^le-^/(^^)(|j| + b|)! 

The corresponding sum in IIX26II is bounded by 

^==0 |il<^+3fc 

X E X2{x,t)(t)j{A{t),B{t),h,a{t),r]{t),x] 

\p\<l+l-k 

l-l ( \ V2 

< E h'/' E \Mm E F'ip) 

k=0 bl<^+l-fc V \j\<J+3k 



DPV{a{t)) {x-a{t))P 



pi 



1/2 



< 



A;=0 

e 



J\k\ 

I- 

E 

fc=0 



\p\<l+l-k \\3\<J+2,k 



(t)2/3(J + 3A;)! 



E ;,H/Vj^V(J + 3A:+H)!, 



(9.13) 



|p|<«+i-fc 



where /? is independent of J. The last sum in this expression is bounded by 



l+l-k 



l+l+2k 



-3k 



(9.14) 



r=0 



s=3k 



Since s<l + l + 2k<3l + lc^g/h and (7 is small, we have 



where a{g) = flo^/g is smaller than one. Furthermore, the sum (|9.14|) is bounded by 

1 - a{g) 



From this we deduce that ( |9.13D is dominated by a constant times 



„-fe/(47i) l-l t,k/2 fk 



fe=0 



A;! 



V- «(^) 

to h'k\ 



< 2 



< 



-b/(An) 



-b/(8h) 



(9.15) 



provided g is small enough, since a{g) ~ g^^"^ as g ^ 0, which is exponentially small and 
integrable in time. 
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Finally, gathering estimates ( |9.8| ), ( |9.12| ) and ( |9.15[ ), we get the existence of positive 7, 
H , G, and C, such that g < G, l{h) = g/h, and h < H imply 



dtUi{x,t)\\/h < Ce-''/^ 



Proof of Lemma |9.2| : Let f E S he the function that is given in spherical coordinates by 



fix) 



For integers g > 0, the maximum absolute value of this function is 



g/2 | y^^^^^ 



(9.16) 



The function / is a normalized eigenfunction of —A + with eigenvalue E = 2q + d. 
Its norm in the Sobolev space Tig for s > satisfies 



< 



Clis) (11/11 + ||(-A)^/V||) < C^is) (1 + i2q + dr/') 



(9.17) 



for some Ci{s). 

If s > d/2, then (1 + Ikf)-'/"^ is in ^^(iR^). So, by H51der's inequality. 



This and (^J7^ imply that is bounded by ^3(5) (l + (2g + rf)^/^^. Thus, 



max I Yq^rn(uj) I < C.^fs) 



T(g+i) 



1 + {2q + dy/') g-«/2g«/2 



The lemma follows from this by an application of Stirling's formula. 



10 More General Coherent States 



In this section we extend all the previous theorems of the paper to allow initial conditions 
that are certain infinite linear combinations of the 6^. 



Proof of Theorem |3.6 : The strategy is quite simple. Let G C have expansion cp 

Cj(j)j{E,E, h, 0, 0, x), and let 

tjjo{x,0,h) = {Ah{a,r])(f){x) 
be our initial condition. By construction, 

ilJo{x,o,fi) = J2 Cj(f)j{n,n, h, a, T], x). 
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For J > 0, we define 

^j{x,0,h) = J2 Cj(j)jmn, h, a, ?7, x), 

\j\<J 

and denote the approximation that arises from this initial condition by 

= J2 Cj{t,n)(f)j{A{t),B{t),n,a{t),r]{t),x). 
\j\<J{Jfi) 

We then have 

e-**^(^)/^V^o(0,/i) 

Thus, to make the error terms exponentially small in h, we need to consider values of the 
cutoff J that grow to infinity with h in a, suitable way, and we also need to control our 
approximation as a function of J . 

In the proofs of all previous theorems, the dependence of the approximation on / governs 
the estimates on the error terms. The h dependence comes through the different choices of 
/ ~ g/h or / ~ g{T)/h, with T ~ ln(l/^). The set C is chosen to give an exponentially small 
contribution as / ^ oo in the last term of ( |10.1| ) with the choice 

J = ul, (10.2) 

for some u > 0. We need only show that the basic estimates in the proofs above are unaltered 
by the replacement of J hj u g/h, for g small enough. 

We can do this because we have been careful to make the J dependence explicit in all 
the key estimates, such as Corollary 

In the contribution to the error term associated with xi given by ( |6.8|) we adapt the last 
step by using the estimate 

for some constants Co(z^) and Ci(z/). Hence, the remainder of the argument for Lemma |0| is 
the same, with updated constants. Since the constants are modified in a time independent 
way, the long time estimates are also unchanged. 

Consider now the contribution associated with X2 in Lemma |6.3| . We first note that ( |10.2|) 



implies (with a slight abuse of notation) J{1) = J{1) + 3/ — 3 = (i/ + 3)/ — 3 so that we still 
have J{1) g/h. The arguments that rely on the smallness of g to allow us to use of Lemma 
( p. 21 ) remain in force. We thus arrive at ( |6.27D . We deal with it by using ( |10.3| ), exactly as 



above, and obtain exponential decay again in case / = g/h. The long time estimates are also 
valid as the time dependence of the constants is unaltered. This shows that Theorems 



and 3.3 are true with our generalized initial coherent states. 



Theorem |3.2| also holds for these initial states provided we can control the sum in (|7^ 



with J(/) = ul. To do so, we first note that the last two factors of (Of) can be bounded by 
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for some so that they are of order e^^^. This is harmless if g is small enough because 
of the exponentially decreasing prefactor. Next, we use < / — 1 to obtain 

for some constant c(z/). Thus, the sum in ( [7.3[ ) is bounded by 

c'(z/)' ^2 {J + 3kf^'^lk\, (10.4) 

fc=0 

for another constant c'(z/), and where is proportional to CqD^ in (|7.3| ). Since 
we can bound ( |10.4| ) by 



which, again, is harmless for sufficiently small g. 



To prove the validity of Theorem p.4| , insert ( p.2| ) and ( |8.3| ) in the estimates above and 
check that the conclusion still holds. This is straightforward. 



Finally, for Theorem |3.5| to hold, we first must consider the contribution associated to 
Xi^i, which relies on (|9.8|) . Here, ( |10.3|) applies directly. Next, the first contribution from 
X26 is ( 9.11| ). It has the form ( |7.3| ) and yields exponential decay in the same way, for g small 



enough. It remains for us to bound ( |9.13| ). With s < Z + 1 + 2A; < 3/, we use the estimate, 

,/{J^. < (z./ + s)-'/2 (z// + s)^/2 < {(u + 3)iy'^mi^ + 3)g/hy/^ 
in ( |9.14| ). The first factor when multiplied by e^"^/ is of order e^^^^ where c is independent 



of g. The final factor allows us to repeat the argument that led to ( |9.15|) . Hence, for g small 



enough, we get an exponentially small contribution in h and the result follows. I 
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